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Abstract 

In the present paper influence of tensor electric and magnetic polarizabilities on spin evolution 
in the resonance deuteron EDM experiment is considered in details. 

It is shown that besides EDM the electric and magnetic polarizabilities also contribute to the 
vertical spin component P3. Moreover, the electric polarizability contributes to the P3 component 
even when the deuteron EDM is supposed to be zero and thereby the electric polarizability can 
imitate the EDM contribution. It is shown that unlike the vertical component of the spin P3 the 
component P33 of polarization tensor does not contain contribution from the electric polarizability, 
whereas contribution from the magnetic polarizability reveals only when the deuteron EDM differs 
from zero. 

Moreover, it is also shown that when the angle z? between the spin direction and the vertical 
axis meets the condition sin$ = (cos$ = y^), the initial value of P33 appears P33(0) = 0. As 
a result, EDM contribution to the measured signal linearly growth in time starting from zero that 
is important for measurements. 

Therefore, measurement of the P33 component of deuteron tensor polarization seems to be of 
particular interest, especially because the nonzero component P33 appearance on its own indicates 
the EDM presence (in contrast to the P3 component, which appearance can be aroused by the 
tensor electric polarizability, rather than EDM). 

1 INTRODUCTION 

At present time the possibility to measure the electric dipole moment of a deuteron moving in a storage 
ring is actively discussed p]l2]- According to PJE] two types of experiments are discussing now: a), 
deuteron energy is chosen to zeroize {g — 2) precession (g is the gyromagnetic ratio) b). deuteron beam 
velocity is modulated with the frequency Qf close to the (g — 2) precession frequency Q, as a result, 
this makes possible to observe the EDM signal as growth in time of the vertical spin component P3I2]- 

These methods can provide for EDM measurement the sensitivity ~ 10~ 29 e • cm. Theoretical 
description of the experiment p] [2j is being done on the base of the Bargman-Myshel-Telegdy equa- 
tion. But as it is shown in [31-0 the Bargman-Myshel-Telegdy equation can not describes deuteron 
spin behavior in such experiments. It turns out that the BMT equation for a deuteron should be 
supplemented with several additions, which describe interaction of deuteron electric and magnetic 
polarizabilities with the electric field in the storage ring and deuteron birefringence in matter. 

Owing to the above in the experiments planned for the deuteron EDM search the contributions 
aroused by the tensor electric and magnetic polarizabilities of deuteron as well as the spin-dependent 
amplitude of forward scattering by the nuclei of a solid or gas target will be also measured. These con- 
tributions could distort the EDM signal and even bring to wrong conclusions about EDM observation. 
They are the systematic errors for the EDM search, which should be eliminated. 



At the same time the above effects being measured in the experiments for EDM search could be 
even used for more reliable limits assignment for deuteron EDM. 

In the present paper influence of tensor electric and magnetic polarizabilities on spin evolution in 
the resonance deuteron EDM experiment is considered in details. 

It is shown that besides EDM the electric and magnetic polarizabilities also contribute to the 
vertical spin component P3. Moreover, the electric polar izability contributes to the P3 component 
even when the deuteron EDM is supposed to be zero and thereby the electric polarizability can 
imitate the EDM contribution. It is shown that unlike the vertical component of the spin P3 the 
component P33 of polarization tensor does not contain contribution from the electric polarizability, 
whereas contribution from the magnetic polarizability reveals only when the deuteron EDM differs 
from zero. 

Moreover, it is also shown that when the angle 1? between the spin direction and the vertical axis 
meets the condition sint? = y| (cost? = yJ~^), the initial value of P33 appears ^33(0) = 0. As a result, 
EDM contribution to the measured signal linearly growth in time starting from zero that is important 
for measurements. Therefore, measurement of the P33 component of deuteron tensor polarization 
seems to be of particular interest, especially because the nonzero component P33 appearance on its 
own indicates the EDM presence (in contrast to the -P3 component, which appearance can be aroused 
by the tensor electric polarizability, rather than EDM). 

2 Interactions contributing to the spin motion of a particle in a 
storage ring 

As it is shown in [HI-IE] considering evolution of the spin of a particle in a storage ring, when measure 
EDM, one should take into account several interactions: 

1. interactions of the magnetic and electric dipole moments with an electromagnetic field; 

2. interaction of the particle with the electric field due to the tensor electric polarizability; 

3. interaction of the particle with the magnetic field due to the tensor magnetic polarizability; 

4. interaction of the particle with the pseudoelectric nuclear field of matter. 

The equation for the particle spin wavefunction considering all these interactions is as follows: 

ih^- = (h + Vedm + v s + v a + VE ud ) *(t) (1) 

where *&(t) is the particle spin wavefunction, 

Hq is the Hamiltonian describing the spin behavior caused by interaction of the magnetic moment 
with the electromagnetic field (equation Q with the only Hq summand converts to the Bargman- 
Myshel-Telegdy equation) , 

Vedm describes interaction of the particle EDM d with the electric field, 

Vedm = -d(p*B + E)§, (2) 

(3 = ^, v is the particle velocity, c is the speed of light. 

Vg, describes interaction of the particle with the electric field due to the tensor electric polarizability: 

Vg = --&ik(E e ff)i(E e ff)k, (3) 

where is the electric polarizability tensor of the particle , E e ff = {E + (3 x B) is the effective 
electric field; the expression (JHJ) can be rewritten as follows: 

V E = a sEef f ~ a T E 2 eff (Sn E ) , n E = (4) 
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where as is the scalar electric polarizability and ay is the tensor electric polarizability of the particle. 

A deuteron also has the magnetic polarizability which is described by the magnetic polarizability 
tensor Interaction of the particle with the magnetic field due to the tensor magnetic polarizability 
is as follows: 



V, 



B 



-^Pik{B e ff)i(B e ff) k , 



(5) 



where (B e ff)i are the components of the effective magnetic field -B e // = (B — [5 x E); © could 
be expressed as: 



V, 



PtB 2 



eff 



Sn B ) , n B 



B- f3x E 



(6) 



\B-(3x E\ 

where f3s is the scalar magnetic polarizability and (3t is the tensor magnetic polarizability of the 
particle. 

ynuci de-scribes the effective potential energy of particle interaction with the pseudoelectric field 
of the target. 



3 The equations describing the spin evolution of a particle in a 
storage ring 

Let us consider particles moving in a storage ring with low pressure of residual gas (10 -10 Torr) and 
without targets inside the storage ring. In this case we can omit the effects caused by the interaction 

\rnucl 
V E ■ 

Let us consider a deuteron moving in a storage ring. According to the above analysis spin behavior 
of a deuteron can not be described by the Bargman-Myshel-Telegdy equation. The equations for 
particle spin motion in condition when the fields E and B are orthogonal to the particle velocity v 
were obtained in They can be written as follows: 



dP 
dt 
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T ft 
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where m is the mass of the particle, e is its charge, P is the spin polarization vector, P{k is the spin 
polarization tensor, P xx + P yy + P zz = 0, 7 is the Lorentz- factor, j3 = v/c, v is the particle velocity, 
a = {g — 2)/2, g is the gyromagnetic ratio, E and B are the electric and magnetic fields in the point 



of particle location, E, 



eff 



(E + fix B), B, 



eff 



(B - (3 x E), n E 



E+pxB 

\E+(3xB\ ' 



n B 



B-/3xE 
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PiknE,k, n ' B i = Pu^Bh Q r (d) are the components of the vector £l(d) (r = 1, 2, 3 correspond to 



x,y,z, respectively). 



n(d) = n + n d , 



fid- d 
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The equations for particle spin motion Q can be rewritten as follows: 

dP 



dt 



[P x £l(d)\ + Q T [n E x n' E ] + Q%[n B x n' B ], 



dP 



ik 



dt 



-(e jkr Pijtt r (d) + e jir P kj Q r (d)) + Vt' T {[n E x P]in Ek + n Ei [n E x P] fc ) + 

+Qy x P]in B k + n Bi [n B X P] fc ) 



(8) 



where 



T 
M 



2 "T^jf Q/ _ 3 ggjw/ q/ 

3 ft ' T ~~ 2 ft ' T 



2 Pt B Ijj_ qIh 



3 PrBjfl O'^ 

T 3 ft ' "T ~~ 2 ft ' "T ~~ 3 " T - 



Suppose that the external electric field in the storage ring E = and a particle moves along the 
circle orbit. 

Let us now consider the equation (JSJ) in the coordinate system that rotates with the frequency 
of particle velocity rotation. In such a system spin rotates with respect to the momentum with the 
frequency determined by (g — 2). The coordinate system and vectors v,E,B as shown in figure and 
denote the axes by x, y, z (or 1, 2, 3, respectively). 
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z 
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Figure 1: 



Therefore, the components of the vectors are: 
P=(P 1 ,P 2 ,P 3 ), 



n E = (1, 0, 0) , n' Ei = PunEl = Pa 

[n E x n' E } 1 = 0, [n E x n' E } 2 = -P31, [H E x n' E ] 3 = P 2 , 

[p x n] 1 = qp 2 , [P x n] 2 = -opi, [P x n} 3 = p 2 , 

[H E x P] x = 0, [H E x P] 2 = -P 3 , [n £ x P] 3 = P 2 , 

n = ^5 = (0,0,0), 

n s = (0, 0, 1) , n' Bi = PuriBi = P i3 

[H B x n' B } 1 = -P 23 , [H B x H' B ] 2 = -P13, [n B x H' B ] 3 = 0, 
[n B x P] x = -P 2 , [n B x P] 2 = P l5 [fi B x P] 3 = 0. 



0) 



(10) 
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Substituting to the system Q we obtain: 

^ = np 2 - n»p 23 , 

*g = -np 1 + - n T )p 13 + u d p 3 , (ii) 

*g = Q T P 12 - cu d P 2 , 



where uj d = — 4 — . 



^ = 2nP l2 + 2u d P 2 3, 



dt 

^ = _ 2 ftP 12 , (12) 



(p u - p 22 ) - n' T p 3 + ^Pis, 
^ = op 23 + n' T p 2 - n'£p 2 - uj d p 12 , (13) 



— np 13 + n'^p 1 - oo d {p 22 - p 33 ). 



^ = -op 13 + njf 

Remember that Pu + P 22 + P 33 = and Pu- = P^i- 



4 Contribution from the EDM and tensor polarizabilities to 
deuteron spin oscillation 

Let us consider the system (|11H13|) more attentively. 

Suppose that deuteron has neither EDM no tensor electric polarizability: in this case the system 
(jl 1H13|) can be expressed: 

* = 
^ = o 

^ = 2^12, 

^ = _20P 12 , (15) 
dt — u > 
^ = -fi(P u -P 22 ), 

^ = ^23, (16) 
^ = -^13- 



This is the conventional system of BMT equations that describes particle spin rotation with the 
frequency equal to (g — 2) precession frequency f2 = ^-B. The component P 3 of vector polarization 
in this conditions is equal to constant = 0) along with the component P33 of tensor polarization 
(^=0). 

Suppose the deuteron EDM differs from zero. Then the above system of equations converts to: 

dEi = -nP 1 +uj d P 3 , (17) 



dP; 
dt 



= 2nP 12 + 2to d P 23 , 
^ = -2QP 12 , (18) 

* = -2^3, 



^ = -n(p u -p 22 ) + UJd p 13: 



dt 

^ = flP 23 - oj d P 12 , (19) 

^ = -nP 13 -u; d (P 22 -P 33 ). 

From (|17B19p it follows that presence of the nonzero EDM makes the vertical component P3 of vector 
polarization oscillating with the frequency of {g — 2) precession Q. 

According to the idea P these oscillations can be eliminated if the deuteron velocity is modulated 
with the frequency Q: 

v = vq + Sv sin (Qt + (pf) (20) 

here ipt is the phase of forced oscillations of the velocity. 
As E e ff depends on f3 = v/c it also appears modulated: 



E eff = E° eff + SE eff sin (fit + ip f ) (21) 

Therefore uj d = E f? f is also modulated with the same frequency. This makes the product 
u d P 2 ~ sm 2 (Qt + iff). Therefore averaging this value over the period of (g — 2) precession gives 
the result time-independent (i.e. = const) and P 3 (t) = P 3 (0) + const ■ t. For better measurement 
conditions it is important to make P 3 (0) = 0. This is the reason to chose particle spin to be in the 
horizontal plane. 

All the above reasoning makes ~ const, too. Therefore, -P33 also linearly grows with time 
^33^) = ^33(0) + const ■ t. However, if the spin lays in the horizontal plane P 33 (0) 7^ 0. 

It is important to note (see below the section I4.2|) that if the spin orientation corresponds to 

cos 2 i9 = I (cosi? = y^> sint? = y^|) then the component ^33(0) = 0, while P 3 (0) / 0. 

Therefore taking $ corresponding to cosi? = \J~^ we can use the component P33 for EDM mea- 
surements, too. 

Let us consider the contribution from the electric and magnetic tensor polarizabilities. Then 
instead the system ()17II19|) we should consider the system (|11I13J1 
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*g = -np 1 + - n T )p 13 + ^p 3 , (22) 
*g = n T p 12 - u d p 2 , 

^ = 2QP 12 + 2io d P 23 , 

= -2QP 12 , (23) 

^ = -n (p n - p 22 ) - n' T p 3 + to d p 13 , 

^ = op 23 + n' T p 2 - n'»p 2 - uj d p l2 , (24) 

4fr = -QP 13 + ^P± - u d (P 22 - P 33 ). 

Some interesting implications follow from I|22l24j) . As it was already mentioned above in the experi- 
ments for EDM search it is planned to measure growth of the vertical component of the polarization 
vector P 3 . 

According to (|22|) time dependence of the vertical component of the vector polarization P3 is 
described by the equation 

dPt 

-± = Q T P 12 - uj d P 2 (25) 

As it can be seen the time dependence of P 3 is determined by both the EDM and tensor polar izability 
of deuteron. It is interesting that the derivative of the tensor polarization component P33 does not 
contain contributions from tensor electric polarizability and is proportional to the EDM only: 

^ = -2u, d P 23 (26) 
at 

Therefore, it is important to measure the component P33, too. According to the above spin orientation 
for this case is determined by the condition cos 2 •& = 3 . 

4.1 Contribution from the tensor magnetic polarizability to deuteron spin oscil- 
lation 

Contributions to spin rotation and oscillations from EDM and polarizabilities are small. Therefore, 
they, being analyzed, could be considered as perturbations to the full system (I22ll24fi and the role of 
each could be studied separately. 

The system of equations considering contribution from the tensor magnetic polarizability [3t is as 
follows: 

4EL = np 2 -n»p 23 , 
dp>=-np 1 + n^p 13 , 

(27) 

^L = nPo 3 -n'^p 2 , 



^ = -np 13 + n%p x 
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Introducing new variables P + = P\+ iP 2 and G+ = P13 + ^23 arid recomposing equations (J2*T|) to 
determine P + and G + we obtain: 

^ = -«QP + + ^G + , 



dt 

or 



dt ~ T ' 

Let us search P + , G + ~ e laJ * then ((25)) transforms as follows: 



p_j_ — r2p_[_ — 



The solution of this system can be easily find: 



T T 

that finally gives 



(u - n) 2 - n£n% = o (28) 



wi, 2 = fi ± y^r^r ( 29 ) 

Rewriting the solution 

P + (i) = c ie - iuJlt + c 2 e-^ 2 * = \ Cl \e- i{uJlt - 8l) + \c 2 \e- i( - W2t - 52) (30) 

Therefore, 

Pi(t) = | ci | cos(a;it — 5i) + | C2 1 cos(u;2i — ^2) (31) 

P2W = — |ci| sin(wit - 61) - \c 2 \ sin^i - 5 2 ) (32) 

According to (|31l32j) the nonzero deuteron tensor magnetic polarizability makes the spin rotating 
with two frequencies u\ and uj 2 instead of O and, therefore, experiences beating with the frequency 

Au; = lu 1 -u; 2 = 2^9^ = ^s££. 

Let us recall now that EDM interaction with the electric field makes the deuteron spin rotating 
around the direction of this field and leads to appearance of P3 component proportional to P 2 (t) 

^ ~ -uj d P 2 (33) 

Therefore, 

dPz 

— — = uj d (|ci| sin(wit - <5i) + |c 2 | sin^^ - fa)) (34) 
dt 

According to the idea to measure the EDM the particle velocity {v = vq + 5vsva.{Qft + (ff)) 
should be modulated with the frequency Qf close to the frequency Q of {g — 2) precession. 
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If the magnetic polarizability is equal to zero, then oj\ = UJ2 = and spin rotates in the horizontal 
plane with the frequency In this case velocity modulation with the same frequency £lf = fl, gives 



^ ~ sm\nt) (35) 
dt 



and the vertical component P3 linearly grows with time. 

However, uj\ 7^ uj% and velocity modulation, for example, with the frequency Q, = u\ provides for 
slow spin oscillation with the frequency uj\ — UJ2 instead of linear growth. 

According to the evaluation jS] the tensor magnetic polarizability (3? ~ 2 • 1CP 40 , therefore the 
beating frequency Aw ~ in the field B ~ 10 4 gauss. 

Measurement of the frequency of this beating makes possible to measure the tensor magnetic 
polarizability of the deuteron (nuclei). 

Thus, due to the presence of tensor magnetic polarizability the the horizontal component of spin 
rotates around B with two frequencies u%, 102 instead of expected rotation with the frequency S7. 

This is the reason for the component P3 caused by the EDM to experience the similar oscillations. 
Therefore, particle velocity modulation with the frequency O provides for eliminating oscillation with 
frequency, but P3 oscillations with the frequency Au> rest (similarly P33). Study of these oscillations 
is necessary because they can distort the EDM measurements. 

4.2 Contribution from the tensor electric polarizability to deuteron spin oscilla- 
tion 



Let us consider now contribution caused by the tensor electric polarizability. From the system (|12j) it 
follows 

(36) 

= _ n d(Pn-P 22 ) _ n>T dF± = _ {AQ 2 + V T tt T )P 12 . 

Thus we have the equation 

^+o,? 2 P 12 = (37) 

where W\% = y^fi 2 + Qt^'t ~ 20, because Qt^'t ^ 

The solution for this equation can be found in the form: 

P12 = ci cos uiut + C2 sin uut (38) 

Let us find coefficients c\ and c 2 : when t = the equation ()38|) gives c\ = -Pi2(0). The coefficient 
C2 can be found from 

d(Pi2) 



dt 



-(t 0) = LO l2 C2, (39) 



therefore 



From the equation (|13|) 



that 



c 2 = -^(t-0), (40) 
W12 at 



HP 

-> 0) = -U (P u (t - 0) - P 22 (t - 0)) , (41) 
C2 = 7, , (42) 
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and 

P\2 = P\2 (0) cos io\2t sin u>i 2 t (43) 

As a result we can write the following equation for the vertical component of the spin P3: 

^1 = q tPi2 ( t ) = n T [P 12 (0) cos 2flt - Pll (°) ~ P22 (°) sin 2 nt] (44) 

As it can be seen the vertical component of the spin oscillates with the frequency 2Q. 

But it should be reminded that according to the equations (J7J) interaction of the EDM with an 
electric field causes oscillations of the vertical component of the spin with the frequency O. According 
to the idea these oscillations can be eliminated if the deuteron velocity is modulated with the 
frequency £lf that should be taken as close to the frequency O as possible: 

v = vq + 5v sin {tlft + iff) (45) 

As a result E e ff depends on /3 = v/c it also appears modulated: 

E eff = E° eff + SE eff sin (tlft + iff). (46) 

The equation describing contribution from the tensor electric polar izability and EDM to P3 looks 
like OH 

^ = n T P l2 - iv d P 2 . (47) 
at 

As P 2 oscillates with the frequency f2, then the product w^Pi contains the non-oscillating terms and 
contribution to P3 caused by EDM linearly grows with time, if £lf = £1. HVlf ^ Q, contribution to P3 
caused by EDM slowly oscillates with the frequency Q.f — Q.. 

It is important that modulation of the velocity v = v + Sv sin (£lft + iff) results in oscillation of 
E 2 e ^ also oscillates with time and appears proportional to sin 2 (fljt + iff). As a result the contribution 
to -P3 caused by the tensor electric polarizability can be expressed as: 

^ ~ AO T sin 2 (^ / t + V 9 / )[P 12 (0)cos2m- ~ P22 ^ sin20f] (48) 

i.e. 

^ ~ - -AO T cos (2n f t + 2cp f )[P 12 (0) cos 2Qt - Pl1 ^ ~ Pa2 ^ sin2nt] (49) 

According to for a partially polarized deuteron beam the derivative depends on the 
deuteron polarization components P\ 2 and Pll (°) P22 

For simplicity let us consider a deuteron beam in pure polarization state. In this case the com- 
ponents P\ 2 and ■ Pll ( )~- F>22 can De written using the explicit expression for the spin wavefunctions. 
Suppose ip) is the unit vector directed along the deuteron spin (1? and (p are the polar and azimuth 
angles (see FigU)). So the spin wavefunction that describes the deuteron spin state with the magnetic 
quantum number m = 1 can be expressed as follows (in the Cartesian basis): 

a x \ I cos 1? cos (p — i sin ip 

Xi ($,</?)= I cL y = — cos # sin if + i cos ip I (50) 



V2 



siwd 



Polarization vector can be written as 

P = (S) = xfSxi = i[a x a*} (51) 
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and components of polarization tensor 



3 2 

(Pik) = xfQikXi = ~^{ a i a l + a k<4 - g}, (52) 



where Qifc is the spin tensor of rank two. Therefore, 

P12 = ^ sin 2^9 sin 2 ??, (53) 
P11-P22 = 3 cos2ygin 2^ (54) 

P 33 = -^(sin 2 ^-0. (55) 

Using (49,53,54) one can obtain: 
dP 3 

r Aft T sin 2 #cos(2fyi + 2<p#) x [sin 2ip cos 20i - cos 2tp sin 2J2i] (56) 

dt 8 

From it follows that ^& slowly oscillates with the frequency (0/ — fi) 

In the ideal case, when Qf = $7 (as it is proposed in for EDM measurement) Q56|) converts to 

^ = -|aO t sin 2 t?cos(2fit + 2p/)sin(2fit-2<p) (57) 
dt 8 

In the general case, when the phases iff and 92 are arbitrary, (|57j) contains terms that do not depend 
on time and, therefore, P3 linearly grows with time like the signal from the EDM does. 
It is interesting that making ipf = — ip gives 

dP 1 

~ cos(2fit - 2(p) sin(20t - 2<p) = - sin(40i - A(p) (58) 

that makes this contribution to P3 quickly oscillating and depressed. But even in this ideal case it rests 
the contribution caused by the tensor magnetic polar izability (JSJ), in real situation 0, ^ Qf, though. 

Measurement of these contribution provides to measure the tensor electric polarizability. 

According to the evaluations j6, olt ~ 10 -40 cm 3 , therefore for the field E e tt = B ~ 10 4 gauss 
the frequency Qt ~ 10 -5 sec -1 . When considering modulation we should estimate AO^ ~ Qt{^) 2 , 

then suppose (^) ~ 10~ 2 - 10~ 3 we obtain Aft T ~ 10~ 7 - 1(T 8 sec , that exceeds the magnitude 
of LOd for the deuteron EDM d = 10~ 29 e • cm. 



5 Conclusion 

In the present paper influence of tensor electric and magnetic polarizabilities on spin evolution in the 
resonance deuteron EDM experiment is considered in details. 

It is shown that besides EDM the electric and magnetic polarizabilities also contribute to the 
vertical spin component P3. Moreover, the electric polarizability contributes to the P3 component 
even when the deuteron EDM is supposed to be zero and thereby the electric polarizability can 
imitate the EDM contribution. It is shown that unlike the vertical component of the spin P3 the 
component P33 of polarization tensor does not contain contribution from the electric polarizability, 
whereas contribution from the magnetic polarizability reveals only when the deuteron EDM differs 
from zero. 

Moreover, it is also shown that when the angle i9 between the spin direction and the vertical axis 
meets the condition sin# = w| (cost? = -wg), the initial value of P33 appears P33(0) = 0. As a result, 
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EDM contribution to the measured signal linearly growth in time starting from zero that is important 
for measurements. 

Therefore, measurement of the -P33 component of deuteron tensor polarization seems to be of 
particular interest, especially because the nonzero component P33 appearance on its own indicates 
the EDM presence (in contrast to the P3 component, which appearance can be aroused by the tensor 
electric polarizability, rather than EDM). 
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